Density-potential mappings in quantum dynamics 
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In a recent letter [Europhys. Lett. 95, 13001 (2011)] the question of whether the density of a 
time-dependent quantum system determines its external potential was reformulated as a fixed point 
problem. This idea was used to generalize the existence and uniqueness theorems underlying time- 
dependent density functional theory. In this work we extend this proof to allow for more general 
norms and provide a numerical implementation of the fixed-point iteration scheme. We focus on the 
one-dimensional case as it allows for a more in-depth analysis using singular Sturm-Liouville theory 
and at the same time provides an easy visualization of the numerical applications in space and time. 
We give an explicit relation between the boundary conditions on the density and the convergence 
properties of the fixed-point procedure via the spectral properties of the associated Sturm-Liouville 
operator. We show precisely under which conditions discrete and continuous spectra arise and give 
explicit examples. These conditions are then used to show that in the most physically relevant cases 
the fixed point procedure converges. This is further demonstrated with an example. 

PACS numbers: 31.15.ee, 71. 15. Mb, 31.10.+Z 
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I. INTRODUCTION 

The essence of the many-body problem lies in our inca- 
pability of handling the huge number of degrees of free- 
dom of many-particle systems and consequently in our 
inability to determine the many-body states. This prob- 
lem spawned a lot of interest into the question whether 
one can devise a closed set of equations for reduced quan- 
tities which do not involve the explicit solution of the 
Schrodingcr equation and in which the many-body cor- 
relations can be approximated efficiently Pursuits in this 
direction have led to various approaches such as many- 
body Green's function theory [1], density matrix theory 
[3, y| and density- functional theory [J, Q- These ap- 
proaches differ in the complexity of the reduced quan- 
tity which is used to calculate the various observables of 
interest. In this work we will focus on the simplest of 
these variables, namely the one-particle density, and ask 
the question to what extent this quantity determines the 
many-body states. 

Within the framework of time-dependent density func- 
tional theory (TDDFT) [||-|8[ this question is asked for 
the special case that the density operator is linearly cou- 
pled to a scalar potential in the Hamiltonian. This linear 
coupling suggests the possibility of a one-to-one relation 
between the scalar potential and the density and hence 
between densities and wave functions. This fact was in- 
deed proven by Runge and Gross pj for the case that 
the potential has a Taylor expansion in time and with 
the spatial boundary condition that the potential van- 
ishes at infinity [9(. Another issue is whether a given 
density can be produced by some scalar potential. This 
existence question, which is usually referred to as the 
w-representability problem is a more difficult one. The 
existence question is nevertheless an important one since 
it allows the construction of an effective nonintcracting 



system having the same density as the one of an interact- 
ing system and thereby convert the interacting problem 
into an effective nonintcracting one. This procedure is 
known as the Kohn-Sham method and forms the basis of 
virtually all applications of TDDFT. The existence can 
be established under the condition that densities and po- 
tentials are Taylor expandable in time [lfj. This condi- 
tion is sometimes too restrictive as has been discussed in, 
e -g- uM- There are, however, indications that both the 
uniqueness and the existence theorems of TDDFT are 
valid under more general conditions that do not require 
Taylor-expandability. As a matter of fact, we know that 
Taylor-expandability is not a necessary condition for the 
validity of these theorems. A first extension was given 
in [l2| to the set of Laplace-transformable potentials un- 
der the assumption of a groundstate as an initial wave 
function. In [13l | a proof of the Runge-Gross theorem for 
dipolc fields without restriction on the temporal form was 
presented. Recently, Tokatly in [14j has given a rigorous 
proof for an arbitrary potential on a lattice. These find- 
ings demonstrate that the restriction to analytic poten- 
tials in time is not fundamental and we can extend the set 
of p otentials beyond Taylor-expandable ones. Recently 
[15l | we have introduced a new proof of the two basic the- 
orems of TDDFT, i.e. the Runge-Gross theorem [7| and 
time-dependent v-representability theorem [lCJ . We have 
reformulated the question whether a one-particle density 
is uniquely defined by an external potential in terms of 
a fixed-point problem. In this way we were able to lift 
the usual restriction of Taylor-expandable potentials and 
densities. 

Here we extend this proof to allow for norms on more 
general function spaces and provide additional mathe- 
matical details. We focus on the one-dimensional case 
as it allows to use established mathematical methods 
from singular Sturm-Liouville theory [ly]. The Sturm- 



Liouvillc operator associated with the density can be 
classified according to the boundary properties of the 
density in which each class gives rise to specific spectral 
properties. We give an explicit relation between these 
spectral properties and the convergence properties of the 
fixed-point procedure. We show precisely under which 
conditions discrete and continuous spectra arise and give 
explicit examples. We finally provide a numerical imple- 
mentation of the fixed-point iteration scheme for the case 
of periodic densities. 

The paper is organized as follows: In Sec. |TT] we in- 
troduce the density-potential mapping and formulate the 
basic questions of this many-body theory as a fixed-point 
problem. We draw attention to the fundamental inequal- 
ity that will give us the opportunity to derive uniqueness 
and existence of a fixed point. In Sec. Mil we will de- 
rive in a general fashion the first part of the afore intro- 
duced inequality by using linear response theory. Then 
in Sec. [IV] we will deduce the second part of the basic 
inequality by using Sturm-Liouville theory. With this we 
show in Sec. IVl uniqueness and existence of a fixed point. 
In Sec lVIl we focus on periodic densities, derive the ex- 
plicit form of the fixed-point iteration and show that an 
elementary numerical implementation of the proposed it- 
eration converges. Finally we conclude in Sec. I VIII 



II. DENSITY-POTENTIAL MAPPING AS A 
FIXED-POINT PROBLEM 

In this section we will introduce the fixed-point for- 
mulation of the density-potential mapping. We will for- 
mulate everything for simplicity in the one-dimensional 
case. Note, however, that the reasoning is independent of 
the dimension of the space in which the particles move. 
All considerations carry over to higher dimensional cases 
unchanged. 

The basic equation we want to examine is the non- 
relativistic equation of motion for a given initial state 
|\&(to)) = l^o) of N interacting particles, i.e. the time- 
dependent many-body Schrodinger equation (TDSE), 

i^|*(t)) = £(M,t)|*(t)>. (l) 

The Hamiltonian in atomic units (e = K — m = 1) is 
given by 

H([v],t)=f + V([v],t) + W, 

where the kinetic energy operator reads 

the interaction energy operator is 

W = -^2 [[dxdx , w(x-x')ft(xcr)ft(x'(T') 

2 a,c' JJ 

x tj}(x'a')-il}(xa) 



and 



V([v];t) = / dxh(x)v(xt) 



(2) 



is the external energy operator with the density opera- 
tor h(x) = ^2 a ijj'(xa)'i(}(xa). The operators if}'(xa) and 
tp(xa) are the usual creation and annihilation field op- 
erators for the spin a and w(x — x') is the interaction 
potential. 

Keeping the initial state |^o) fixed for all further con- 
siderations, we observe that there is a mapping between 
external potentials v(xt) and the time-dependent wave 
functions, i.e. each external potential generates an asso- 
ciated |^([w],i)) by propagation of the Schrodinger equa- 
tion. We assume that for every v <S 03, where 03 is the 
set of potentials under consideration, the Schrodinger 
equation (TTJ) has a unique square-integrablc solution 
| "P ([«],£)). Actually, the wave functions are at least spa- 
tially two-times (weakly) diffcrentiable (171 ) . 

Since the wave functions are uniquely defined by the 
potential we immediately find that also all expectation 
values are uniquely determined by the potential, i.e. for 
a physical observable represented by an operator O we 
have 

o(M,*) = (*([«], t)|6|*(M,t)). 

The density is a special observable since it couples di- 
rectly to the scalar potential in the Hamiltonian as is 
directly clear from Eq. ^. Hence we may expect a one- 
to-one relation between densities n([v], xt) and potentials 
v(xt). However, since physical observables are gauge in- 
variant such a one-to-one relation can only be expected 
up to a trivial spatially constant shift c(t) in the poten- 
tial. This would mean that there are no two potentials 
differing more than a gauge that generate the same den- 
sity. If this is true, then the density uniquely determines 
the potential, i.e. v([n],xt), and we find following the 
above reasoning, that the wave function (up to a physi- 
cally irrelevant phase factor) is uniquely determined by 
the density, i.e. |^([n],t)). Consequently all observables 
become functionals of the density and we can in principle 
calculate all quantum mechanical expectation values by 
only knowing the density of the system. 

Our first task now is to determine a way to verify that 
the potential is determined by the density alone. To do 
so, a direct relation between both entities is desirable. 
The obvious way is to use the Schrodinger equation and 
deduce such an interrelation. So we start by the evolu- 
tion of the density, which is controlled by the Heisenberg 
equation of motion. This leads to the well-known conti- 
nuity equation 



d t n(xt) = -d x j(xt), 



(3) 



where d x = d/dx and similarly for the time variable. 
Here j(xt) is the expectation value of the current-density 
operator which is defined by 

j( x ) = 77- ^2 [V' t ( xcr ) \dxi>{xo-)\ - (d x ip\xa)J %j}{xa) 



In order to make the dependence on the external poten- 
tial explicit, we apply the Heisenberg equation on the 
current-density operator and find 

d t j{xt) = -n{xt)d x v{xt) - (d x T xx (xt) + W x {xt)), (4) 

This equation describes the local-force density of the sys- 
tem. The momentum-stress tensor T xx (xt) is defined to 
be the expectation value of the operator [lj| 

f xx (x) = ^{(c^Ora)) dj(xa) 

a 

-\dl^{xa)^{xa))), 

and the divergence of the interaction-stress tensor W x (xt) 
is found as the expectation value of 

W x (x)=^2 fdx'(d x w(x-x'))^(xa)4>Hx'a') 

a, a' 

x i]j(x'a')il>(x(j). 

By using the continuity equation ((3j) in Eq. (|4]) we find an 
explicit relation between the density and the potential: 

-d x [n([v],xt)d x v(xt)] =q([v},xt)~d'?n([v],xt). (5) 

In this equation 

q([v],xt) = (*([v},tMx)\*([v],t)), 

q(x) = d x (d x f xx {x) + W x (x)). 

Equation ([5]) enables us to investigate the density to po- 
tential mapping n <— > v. We can do this by inserting a 
given density into the equation such that we obtain 



d x [n(xt)d x v{xt)} = q([v],xt) — d^n(xt). 



(0) 



We can then search for a potential v that solves this 
equation. To do this we need also to give the initial 
state in order to calculate q([v],xt) on the right hand 
side of the equation. The are now two cases to consider. 
In the first case we assume that n(xt) = n([u],xt), i.e. 
it is the density obtained by propagation of the TDSE 
using some potential u and the given initial state. We 
then know that v = u is a solution to this equation. 
If there is no other potential that solves the equation 
then there is a one-to-one correspondence between the 
density and the potential. The Runge-Gross theorem is 
therefore equivalent to the uniqueness of a solution of 
Eq. ([6]). In the second case, we consider a density of 
which we do not a priori know whether it is obtained 
from a time-propagation of the TDSE. In this case there 
are two possibilities to consider. The first possibility is 
that there is no solution to Eq. ([6]). In that case the given 
density is not v-represcntablc for the given initial state. 
The second possibility is that we find a solution v. In 
that case it is not immediately clear that this potential 



v produces the given density. However, if we propagate 
the TDSE using this potential we satisfy Eq. (|5|) in which 
q{[v], xt) by construction is the same as in Eq. ^j. If we 
therefore subtract both equations we obtain 

dfp(xt) - d x [p(xt)d x v(xt)} = 

for the density difference p(xt) = n([v],xt) — n(xt). For a 
given v this a linear and homogeneous differential equa- 
tion for p(xt). Let us now discuss its initial and boundary 
conditions. It follows immediately from the equation of 
motion of the density operator that the density from the 
time-propagation satisfies the conditions 



n(xt ) 
d t n(xt)\ to 



(*oKa:)|tto>, 
-(#o|3xJ»I*o) 



(7) 



Also the given density n(xt) must satisfy these condi- 
tions otherwise we obviously can not find a potential v 
producing this density and the chosen density would not 
be w-representable. If we therefore choose n{xi) to have 
these initial conditions then p{xt) satisfies 

p(xt ) = d t p{xt ) = 0. 

Furthermore to have a solution the density n(xt) must 
have the same spatial boundary conditions as n([v],xt) 
which are dictated by the TDSE. This gives two further 
conditions. Finally we could add the additional condi- 
tion that p(xt) integrates to zero when integrating over 
space. Since we have already five conditions on a linear 
differential equation that is second order in time and first 
order in space it is clear the solution p(xt) = is the only 
one. We thus find that if v is a solution to Eq. ^ then 
n{xt) = n([v],xt). 

After having fixed the boundary conditions for the po- 
tentials under consideration we can then ask the question 
whether there is any other solution u/d generating the 
same density. This means we want to examine whether 
there is another u G 23 for which Eq. © holds such that 
by subtraction we would find 



d x [n(xt)d x u)(xt)} = q([v],xt) - q([u],xt). 



(8) 



where we defined u>(xt) = v(xt) — u(xt). In the case of 
the original Runge-Gross proof this question is answered 
by taking repeated time derivatives of this equation at 
the initial time to. This assumes that all time deriva- 
tives of n{xt) and to(xt) in to exist. At time to we have 
q([v),xt ) = q([u],xt ) = (^ \q(x)\^f ) and hence 



d x [n(xt )d x uj(xt )} = 



(9) 



Since u and v have the same boundary conditions u{xt) 
vanishes at the boundaries and the unique solution is 
uj(xto) = 0. If we now take the first time derivative of 
([8]) we obtain the equation 

-d x [n(xt )d x d t oj(xt )} = 

dy(^o\[q(x),n{y)]\^o}u(yt ) 



-d x [d t n(xt )d x u}(xt )} 



Since the right hand side vanishes we find for dtu>(xto) 
the same equation Q as for uj(xto) and we find that 
dtUj(xto) = 0. Continuing this way we find that all 
time derivatives of u(xt) vanish at the initial time to, 
i.e. d^ui(xt ) = for all integers k > 0. We thus see 
that it is a necessary condition for two potentials to 
give the same density that all the time derivatives of 
their difference in to vanish. Therefore if one of those 
derivatives for some k does not vanish the two poten- 
tials can not give the same density. However, it is still 
possible that all dtU>(xt ) = while the potentials u 
and v are still different. For example, a function of the 
form uj(xt) = f (x) cxp (—1 / (t — t ) 2 ) has all its time- 
derivatives vanishing in to- To eliminate such cases we 
have to demand that the function u>(xt) is equal to its 
Taylor expansion around to- With this additional condi- 
tion the vanishing of all d^uj(xto) implies that u>(xt) = 0. 
Therefore the mapping from the set of Taylor expandable 
potentials around to with the given boundary conditions 
to the set of densities produced by it is one-to-one. This 
is the statement of the original Runge-Gross theorem. 

Now one can pose the question, whether the Taylor- 
expandability of the external potentials is an essential 
condition for the Runge-Gross theorem to hold. To 
rephrase, can we prove the one-to-one mapping between 
the potentials and densities also for other, possibly more 
general sets 93? As a matter of fact, we know this to be 
true. As discussed in the introduction there are already 
several extensions [12H141 ] which demonstrate that the re- 
striction to analytic potentials in time is not fundamen- 
tal. In this work we want to broaden the set 2J of allowed 
potentials even further. We will do so by using Eq. ([6]) 
in order to define a mapping we will call J 7 , which maps 
potentials to potentials. The first part of this mapping 
concerns the right hand side of Eq. © : we take &v ^%3 
and propagate the initial state |<l>o) in a fixed finite time 
interval [ioj^l with this potential. From the associated 
wave function |\I/([i>o], £)) we calculate q([vo],xt). This 
procedure we denote by 

V : vo H> q[v ]. 

It calculates for every potential the corresponding diver- 
gence of the internal- force density. The second step of the 
mapping T we identify by the left hand side of Eq. ([5]): 
we take the previously determined q([vo],xt) and solve 
the linear differential equation (with the previously cho- 
sen boundary conditions) 

- d x [n{xt)d x vi{xt)\ = q([v ],xt) - d 2 t n{xt), (10) 

in order to calculate a new potential v\. This operation 
we designate by 

V : q[vo] i-4 V\. 

It computes the potential to a given divergence of the 
local force-density and a chosen one-particle density n. 
Next we define the combined map (see Fig. [1} 




T:v ^(VoV)[vo]=v 1 , 



(11) 



FIG. 1: The potential-potential mapping T of Eq. (|11|) as 
composition of the mappings V and V. 



which transforms our original potential vq into V\ . What 
is the connection between T and the original problem of 
Eq. ([6]) or synonymously Eq. ([8])? There are now two 
cases to consider. In the first case we take the density in 
Eq. ^ to be one coming from a potential v, i.e. n = n[v], 
then clearly 

J-[v] = v. 

Hence, T maps solutions of Eq. ([5]) to itself. Thus so- 
lutions of the original equation are fixed points of the 
mapping J- '. Therefore, instead of asking whether there 
is a unique solution to Eq. ([5]) we can equivalently ask if 
T has a unique fixed point. 

In the second case we can insert a density into Eq. © 
for which we do not a priori know whether it can be 
generated by a potential v. Then the existence of a fixed 
point v guarantees, as shown above, that the density is v- 
representable, i.e. n = n[v] . Therefore existence of a fixed 
point is equivalent to the iJ-rcprcscntability of the density. 
The u-representability question is essential for the exis- 
tence of a Kohn-Sham system in density functional the- 
ory, since there we ask whether a given density obtained 
from an interacting system can also be reproduced in a 
nonintcracting system. Therefore the Kohn-Sham sys- 
tem only exists when the w-representability question can 
be answered positively. 

A first trivial test of u-represcntability is that the 
given density n(xt) has to fulfill the initial conditions, 
i.e. Eq. ([7]). Hence, we have to have an appropriate 
initial state with the right density. Only then Eq. (|6]) 
can have a solution v at all. In order to investigate 
existence of a solution of Eq. ^ we again have to fix a 
boundary condition. Assuming Taylor-expandability in 
time of both, the potential as well as of the density, it 
was shown in [l0( how to construct the unique potential. 
This proof of u-represcntability complemented the 
original Runge-Gross proof and lent justification to the 



time-dependent Kohn-Sham scheme. Here we want also 
to go beyond the assumptions of the original extension 
of the Runge- Gross proof and use a formulation in terms 
of a fixed-point problem. 



where v\ = vq + XAv and Av = v\ — vq. A straightfor- 
ward calculation in the interaction picture of quantum 
mechanics and expanding the associated evolution oper- 
ator in powers of e leads to 



In order to answer the raised fixed-point questions 
we will apply the following inequality: 



dO([v x ],xt) 



\\T[vi] - T[v ]\\a < a\\vi - v \\ a 



(12) 



with a < 1 and where || • || Q is an appropriate norm de- 
pending on a positive parameter a on the space of poten- 
tials. This inequality will directly guarantee uniqueness 
of a given fixed point and with some further work we 
can deduce under which conditions a fixed point exists. 
In order to derive inequality (|12[) we will first obtain an 
inequality of the form 



C 
\q[vi] - q[v }\\ a < — f=||ui - i>o| 



(13) 



where C is a positive constant and where p > 1 is a con- 
stant determining the function space. We subsequently 
derive the inequality 



||.FM - J>o]|U < D\\q[ Vl ] - q[v }\\ c 



(14) 



for a positive constant D. Using these two inequalities 
(JT3J) and (fT4")) we can immediately construct the required 
Eq. (fj"2"j) where a = CD/ ^a. If we choose the positive 
parameter a > (CD) P then clearly a < 1. We note here, 
that one could in principle also use different norms for 
the space of potentials and for the g-functions. However, 
for simplicity, we keep those spaces identical. 



III. GENERAL LINEAR-RESPONSE 
INEQUALITY AND THE a-NORM 

The major ideas for inequality (I12[) and the subsequent 
fixed-point approach are found in the derivation of in- 
equality ()13l) . It will not only introduce the afore men- 
tioned a-norm, which is the most important ingredient 
for making the proof work, but due to its universality, the 
derivation allows for different explicit realizations. Thus 
we can formulate the derivation for a general operator 
O and can keep the associated function spaces undeter- 
mined for the moment being. Again, the derivation ap- 
plies directly to the three-dimensional case. The general 
idea that we present is a very simple one. We just want 
to quantify the physically intuitive idea that if two poten- 
tials ^o and vi are close then also the expectation values 
0[v\] and 0[vq] calculated from them by time-evolution 
of the TDSE are close in some norm. We start by calcu- 
lating the non-equilibrium linear response of an operator 
O with respect to a parameter A £ [0, 1], i.e. 



dO([v x },xt) 
dA 



= lim 



0([v x + eAv],xt) - 0([v x ],xt) 



dA 

t 



df / dx'{* \ [Oh x (xt),n Hx (x't')}\^ }Av(x't'). 



where [. , .] is the usual commutator. The subindex H\ 
indicates the operators in the Heisenberg picture for the 
Hamiltonian H(jv\],t), i.e. 

<W*) = U([vx};t 0> t)dU([v x };t,to), 

with U([vx]',t,to) the unitary evolution operator associ- 
ated with H([v\],t). Since the Hamiltonian H([v\],i) 
is explicitly time-dependent this evolution operator is a 
time-ordered exponential. Then by the fundamental the- 
orem of calculus [13 we arrive at 



O([ vi },xt)-O([vo},xt) 



M dO,. . , 

dA— KW 
dA 



dt' / dx'xixt^'t^faix't') - v (x't')), 



where we have defined 



x (xt,x't') = -i / dA(* |[0 Hx (»t),fifl- A (a/* / )]|*o)- 

Jo 

(15) 

The linear response kernel \ i s assumed to be bounded 
in some properly chosen function space with norm 



l/WII 



dr\f(rt)\P 



Up 



(16) 



with p > 1. In the following we will always use this norm, 
unless explicitly stated otherwise. It is now straightfor- 
ward to derive that 

||0(H,u-o(M,uir <<>(*) / t ||A«(*')ll p , (17) 

Jt 
where the constant C(t) is the operator norm defined as 



c*(t)=su P "tawii* , 

s #o J! dt'\\g(t')\\P 



(18) 



This constant has an intuitive interpretation; it simply 
compares the norm of g to that of X9 an d searches for 
its largest possible ratio, i.e. the maximum amplification. 
We observe at this point, that \ = x[ v o, vi] as can be seen 
directly from Eq. (|15[) . Thus, the linear response kernel 
depends on the choice of v\ and vq. As a consequence 
also the operator norm C(t) = C([vo,vi],t) has the same 



dependence. The integral on the right hand side of in- 
equality (|17p can now be manipulated as follows 



dt' \\Av(t')\\ p = / <ft'e- a(t, -* o) e a t*'-*»)||Ai;(t / )|| p 

Jto 

a(t-to) 



"'to 

In this equation we defined the norm [19( 



(19) 



|Ai 



la,* 



sup 
t'e[t ,t] 



-a(t'- 



to) ||A<;(*')ll P 



where a is an arbitrary positive number. Such norms are 
commonly used to prove existence of solutions to differ- 
ential [2(| and integral equations [23, [22| . In Appendix 
fAl we show that all a-norms are equivalent and hence we 
can change a without changing the function space that 
we are considering. We now insert inequality (|19p into 
Eq. fTT]). multiply both sides with e ~ Q( *~ to) and take the 
suprcmum over [to , i] ■ We then obtain 



\0[vi] 



O[vo]\\* < ^\\Av\E t , 



where C(t) — sup t , e r to t i C(t'). If we define 

ll/IU = II/IUt 



(20) 



and C = C(T) we find 



O 



||0[ Ul ]-0[ Vo ]|| a <--|| t ; 1 - Uo || 

on [tg, T\. In this derivation we did not explicitly select a 
function space, i.e. picked a value for p. As we used the 
same norm for the potentials as well as the O-functions 
the choice of p dictates the space of potentials. Thus the 
set of potentials under considerations 23 is part of the 
space of functions v which have finite ||-|| a -norm. If we 
now choose O = q(x) then we have derived inequality 
(|13p within the appropriate a-norm. 



IV. THE V-MAPPING 

In what follows we will derive the second inequality 
and then obtain the main statements of the fixed-point 
approach. In order to get the inequality (JTJ| we will use 
the norm (|16j) with p = 1 for the case of periodic bound- 
ary conditions. For more general cases we use the norm 
with p — 2 since then we can make use of what is called 
Sturm-Liouville theory. However, the inequality may well 
be valid for more general norms and hence for broader 
sets of potentials than treated here. The restriction to 
one-dimensional Sturm-Liouville theory is convenient, as 
the theory and all its strongest statements are usually for- 
mulated for the one-dimensional case only [16J |. For an 
extension to higher dimensions in the context of TDDFT 



the authors provided some results in [23J, |24| . 
The V-mapping is defined by Eq. (|10j) and therefore by 
the solution v\ to the following inhomogeneous equation 
as a functional of its inhomogeneity £ 



where 



d x [n(xt)d x vi(xt)] = C(xt), 



((xt) = q([v ],xt) - 9 t 2 ra(xt). 



(21) 



(22) 



We consider this equation on an interval from a to b. We 
already note that an important property of the function 
C is that it is orthogonal to the constant function, i.e. 



0= / dxC(xt). 

J a 



(23) 



This is a consequence of the fact that q is a divergence 
and that the number of particles is conserved. This fact 
will be important later in our discussion. We can directly 
integrate the Eq. ([21]) and its general solution is given by 

vi(xt)= [ dyG t (x,y)C(yt)+c [ dy-^— + d, (24) 



in which 

Gt(x,y) 



P(l/ 



n(yt) 



e(x-y)} / dz—r, 
J y n{zt) 



where 9 is the Heaviside step function and c and d are 
constants determined by the boundary conditions. Note 
that the last two terms in Eq. (|24l) simply represent 
the most general homogeneous solution (i.e. £ — 0) of 
the differential equation ([2"T]) . We see that in Eq. ([24]) 
the integrals may diverge when the density goes to zero 
at the boundaries. Let us therefore first consider the 
case in which this does not happen. The only physical 
relevant case where this applies is the case of periodic 
systems. In this case periodic boundaries are imposed 
on the Schrodingcr equation and hence on the densities 
and potentials. We can then identify boundary a with 
b and we thus have a finite domain on which we require 
v\(a) = Vi(b) as well as i>i( a ) = v i(fy, where the prime 
means a derivative with respect to the spatial coordinate. 
The solution ([24|) with these boundary conditions is given 
by 



where 



vi(xt) = / dyK t (x,y)((yt), 



r]{xt)r}{yt) 



(25) 



Kt(x,y) = G t (x,y) 



Ja & n(y 



(yt) 



and 



V{Xt) = \ (/ ^y 



(yt) ' J b n(yt) 



Then employing the periodicity of the density we can 
readily check that the kernel K t satisfies the boundary 
conditions 

K t (a,y)-K t (b,y) = 0, 
1 
d x K t (a,y) - d x K t (b,y) 



n(at) 



and similarly in the y- variable. Then as a consequence 
of these boundary conditions and the fact that ( is or- 
thogonal to the constant function (see Eq. (|2"5)l ) we see 
that vi satisfies the required boundary conditions. From 
Eq. ([25)) we now see that 

\vi(xt)\< f dy\K t (x,y)\\((yt)\ 



< max\K t (x,y)\ / dy\C(yt)\. (26) 

y£[a-b\ J a 

Since Kt is a continuous function on a finite domain it 
attains a maximum and we thus see that V\ (xt) is finite 
whenever |£| is integrable. This implies in particular that 
| vi | is integrable itself and that 



fa pb 

dx\ Vl {xt)\ <D t / dy\((yt)\, 



where 

D t = max\K t (x,y)\. 

Now we can use this inequality in our iteration scheme 
and consider the distance between successive potentials 
of Eq. (|5j! in the p = 1 norm of Eq. ((TBI : 

||«2(*)-«i(t)||= / dx\v 2 (xt)- Vl (xt)\ 

J a 

< D t dx\q([vi],xt)-q([vo],xt)\ 

J a 

where we used the explicit form Eq. ((2"2"|) of the inho- 
mogeneity From this inequality it immediately follows 
that 

||«2 - vi\\ a < D\\q[v{\ - q[vo]\\ a 

where we used the a-norm defined in Eq. (|20|) and 

D = max D t . 

te[t a ,T] 

The result above was derived for non- vanishing densities. 
However, if the density vanishes at the boundaries then 
the integral in Eq. (|2"3|) may diverge. To treat this case 
we make use of singular Sturm-Liouville theory since this 
theory allows us to have divergent potentials provided 
that they are square integrable. This then naturally leads 
to the consideration of the case p = 2 in Eq. ((To]) . We 
note that Eq. ((2~Tj) has the form of a Sturm-Liouville 



boundary-value problem (see Appendix [Bj parametri- 
cally depending on t. One can always choose the bound- 
ary conditions in such a way that the Sturm-Liouville 
operator 

S t = -d x [n{xt)d x ] 

is self-adjoint in the Hilbert space of square-integrablc 
functions. Then we can solve Eq. (J10[) and thus properly 
define the mapping V, since we are able to expand the 
inhomogeneity of Eq. (|10p in terms of a time-dependent 
orthonormal eigenbasis {ipi(xt)}, i.e. 

- d x [n(xt)d x cpi(xt)] = \i(t)(pi(xt). 

The eigenvalue Xo(t) = is a special one for which we 
can find the eigenfunction explicitly as 



ipo(xt) = ci(i) +c 2 (t) 



n{yt) 



with c\ (t) and c 2 (t) are constants and xq is an arbitrar- 
ily chosen point a < xq < b. The quantities c\(t) and 
C2(t) are determined by the boundary conditions and the 
normalization. In the following we will always choose 
boundary conditions in such a way that C2(t) = and 
hence tpo(xt) = c\(t) is simply the constant function. The 
appearance of the constant function is a consequence of 
the gauge freedom in the Schrodinger equation, i.e. a con- 
stant shift in the potential will not change the density. 
For a more detailed discussion of the boundary condi- 
tions we refer to Appendix [B] We can now expand C, in 
the cigenfunctions as follows 

oo oo 

C(xt) =X><(st)fo(t)|C(*)) =£C<(*)w(st), 

i=l i=l 

where we used the standard inner product 
(/!<?>= / dxf(x)g(x). 



The trivial zero eigenvalue does not appear in this ex- 
pansion, since the scalar product of £ with the constant 
function is zero as noted before in Eq. (|2"5)l . Therefore 
the solution to Eq. ((TO]) 



vi(xt) = ^jTj7 l Pt(xt) 
i=l '*■ ' 

is perpendicular to the time-dependent constant function 
too. However, we can always add such a constant to 
the unique potential without changing the physics. For 
comparison we note that the solution can be written in 
a form analogous to Eq. (|2"5|) 



v 1 (xt)= / dyT t (x,y)C(yt), 



where we defined the Green's function 12 5H 



oo 1 

T t (x,y) =^jj-r¥>i(xt)(p*(yt). 

i—l l \ ' 

Since Ao = < |Ai| < |Aa| < • • ■ we find that v\ is square 
integrable if C, is, because we have the simple inequality 



n-iWii 2 = E 



at) 



Mi) 



< 



1 



x 1 (ty 



X>(*)i' 



8 = 1 
1 



Ai(t) : 



r||CWH 2 <^. 



Thus the mapping V is well-defined if we assume £(t) 
to be square-integrablc. After we have shown that the 
mapping V is well-defined, we can in a next step use the 
expansion in an eigenbasis to derive the second inequality. 
If we look at two successive potentials, say .Fhvi] = «2 and 
T[vq] = vi, we have in accordance to Eq. © 

-d x [n(rt)d x (v 2 {xt) - vi(xt))] = q([vi],xt) - q([v ],xt). 

We can then expand v-i — v\ and q[vi] — q[vo] in terms of 
the eigenfunctions of St and similarly obtain 



IM*)-t;i(t)||< 



|Ai(i)| 



\q[vi](t)-q[vo](t)l (27) 



If we now multiply Eq. ([27j) with e " Q (*-*°) and take the 
supremum over the interval [to > T] we arrive at inequality 
(flit in which 
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max , 
*e[to,T] l|Ai(i)| 



As pointed out before, in general we find a self-adjoint 
operator St with well-known spectral properties. These 
spectral properties can be related to the behavior of the 
density close to the boundary and a full classification 
of the different cases can be found in the Appendix iDl 
We either will have a discrete spectrum of normalizable 
eigenfunctions or we have a continuous spectrum of gen- 
eralized (non-normalizablc) eigenfunctions or a combina- 
tion of both. In all cases we can expand in those eigen- 
functions since in the case that we have a continuous 
spectrum we can replace the sum over eigenvalues in the 
expansions by an integral (see Appendix |E| . However, 
to derive our inequality we then have to make sure that 
there is a spectral gap between the zero eigenvalue and 
the continuum. We give an explicit example where we de- 
duce a continuous spectrum gapped away from zero and 
the associated generalized eigenfunctions in Appendix [Cl 
In Appendix [D] we show that a spectral gap exists when- 
ever the density does not decrease faster than a quadratic 
function to zero at the boundaries. 

Let us now give an explicit example of a Sturm- 
Liouvillc equation for which we calculate the eigenvalues 
and eigenfunctions in time. The example involves the 




space (a.u 



FIG. 2: (color online). The density evolving in time from 
t = until T = 0.5. 



physical situation in which we have a single particle on 
a ring, i.e. the Schrodinger equation on the interval from 
— 1 to 1 with periodic boundary conditions (the point 1 
is identified with —1). As initial state we choose 



#o(z) = c(exp(-l/(l-x 2 )) + l) 



(28) 



where c is the normalization constant that normalizes 
the wave function to one. We propagate this initial state 
with the external potential 



v(xt) = sin 2 (7ra;) sin(lOt) 



(29) 



in atomic units for a short period of time, say from 
to = to the final time T = 0.5 and calculate the den- 
sity n(xt) and St- We use a Crank-Nicholson scheme 
for time-propagation on an equidistant grid and calcu- 
late the time-dependent density (see Fig. (J3J). By di- 
agonalizing the Sturm-Liouvillc operator with periodic 
boundary conditions we find the eigenvalues as well as 
eigenfunctions in time (see Fig. §5§ and flU). From this 
we deduce the value of the constant in inequality (TH)) 
for the current example as D = max te [ to .T]{|Ai(t)| -1 } = 
l/min tG [ t0jT ]{|Ai(t)|} ~ 0.23. We note that at certain 
times the eigenvalues get degenerate (indicated with an 
arrow in Fig. Q) and therefore at the crossing point the 
cigenfunction with the lowest eigenvalue changes discon- 
tinuously without changing the number of nodes. This is 
clearly visible in Fig. (QJ. For the lowest lying non-trivial 
eigenfunctions (in the periodic case this amounts to two 
nodes) we can see how they change at the indicated times. 
This is a special feature of the periodic case, which can- 
not happen in the case that the boundary condition at 
one boundary does not depend on the boundary condi- 
tion at the other boundary. The calculated eigenvalues 
(see Tab. |T])) and eigenfunctions nicely agree with cal- 
culations done using the SLEIGN2 Sturm-Liouville code 

Before we use the derived inequality (fT2|) to deduce 
uniqueness and existence of a fixed point we will shortly 
summarize our assumptions. For the periodic case where 
the density is non- vanishing we work with the p = 1 norm 



TABLE I: The four lowest eigenvalues of the Sturm-Liouvillc 
eigenvalue problem for to = and T — 0.5 in atomic units 



(i 







T = 0.5 



1st 
2nd 
3rd 
4th 



4.5646 

4.8502 

18.6924 

18.8226 



4.7333 

4.8252 

19.0458 

19.1178 



and hence we require q([v], xt) — dfn(xt) to be integrable. 
This automatically implies via inequality (|26|) that the 
potential is finite everywhere. For the case, of vanishing 
densities we use Sturm-Liouvillc theory and the p = 2 
norm, and we therefore demand q([v], xt) — dfn{xt) to be 
square- integrable in the domain of St, i.e. those v that 
fulfill the chosen boundary and some further regularity 
conditions (they are in D max as defined in Appendix |B| . 
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FIG. 3: (color online). The four lowest eigenvalues in time 
from t — until T = 0.5 for the Sturm-Liouville eigenvalue 
problem. The arrows indicate the time of degeneracy of the 
two lowest lying eigenfunctions. The constant 1/D ~ 4.39 as 
indicated at time t ~ 0.17 




Y'o time ( a - u -) 



FIG. 4: (color online). The lowest eigenfunction in time from 
t — until T = 0.5 for the Sturm-Liouville eigenvalue prob- 
lem. The lowest lying eigenfunction at certain points becomes 
degenerate with the second lying eigenstate, whereafter they 
change position in the spectrum. 



In fact, we can weaken the requirements somewhat, since 
for the inequalities (fl~3|) and (fl"4|) we only need the above 
mentioned properties for the difference potential t>2 — i>i 
as well q[vi] — q[vo], meaning that the potentials at the 
various iterations may have singularities provided that 
they occur at the same spatial points. 



UNIQUENESS AND EXISTENCE OF A 
FIXED POINT 



V. 



In this section we will use inequality (jT2|) to derive 
uniqueness and existence of a fixed-point. The presented 
proofs follow the logics of the Banach fixed-point theo- 
rem [27J . 

First we will show uniqueness of a given fixed-point. 
Let us define by 93 the set of potentials with a given 
set of boundary conditions. Let us take a potential v 
out of this set and calculate n{[v\,xt). This is then, by 
definition, a w-representable density. Let us now assume 
that there is a second fixed-point u in 03, i.e. a potential 
yielding the same density. Then we find by choosing 
tfa = 2 CD for this pair of potentials in Eqs. (|12I14[) 
that 

\\v - u\\ a = \\F[v] - F[u}\\ a < -\\v - u\\ a . 

This can only be true if 



and thus we have u = v. 

Let us now address the existence of a solution to 
Eq. (|6j). This is a w-representability question for a given 
density. Before we present the actual proof we give a 
simplified example that illustrates the physical meaning 
of the assumptions made for the response function \- 
To do so, we look at an analogy in which potentials are 
represented by real numbers v € [a, b] and where an 
observable 0{v) is represented as a real function which 
maps to R. Thus we associate with v a point on the real 
axis. Following the reasoning of Sec. IIIII we then look at 
the derivative of the function O with respect to some 
parameter A e [0,1], where we define v\ = vq + XAv 
with Av = vi — vq. From the fundamental theorem of 
calculus we can thus derive that 



dO[v x ] 



O(wi)-O(«o)= / dA 



n 



dA 



1 dO[v x ]dv x r , 



ii 



di>A dA 
where L d\dO[v\]/dv\ = x[uo,fi]. Then we find that 



C[vq, vi] = max 

VE[Vo,Vl] 



dO\v] 



dv 
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is an upper bound for the slope of the tangent and we 
find that 

\0(vi) - O(v )\ <C[v ,vi]\Av\, 

in correspondence with Eq. (jTTJ) . Thus if V\ approaches 
vq also 0(y\) approaches 0(vq), i.e. 0(v) is continuous 
on [vq, vi] (see Fig. (0). If there is a maximum slope or 
derivative of the function 0(v) when we range over all v 
in [a, b] then there exists a constant 



^'8 



sup C[vo,vi], 



when we range over all pairs i>o and v\ . This means phys- 
ically that small changes in v cannot lead to arbitrarily 
large changes in the observable 0{v). 

We will now show that similar assumptions for the 
response function (|15l) imply the existence of a fixed 
point for T ' . We see from Eq. (fT8|) that the constant 
C = C(T) in Eq. ([15)) is dependent on the response 
function x an d hence via Eq. (| 15[) on potentials vq and 
Vi, i.e. C = C[vq,vi\. We assume that a constant 
C SU p = su P Uo C^c^^o]] exists when we range over all 
potentials vq in the domain 03 of potentials with a certain 
boundary condition. Physically this amounts to assume 
that one cannot induce arbitrarily strong changes in the 
internal forces by weakly perturbing the quantum system 
during a finite time. Following our simplified reasoning 
from above, we essentially presume that the slope of the 
tangent of q[v] does not become infinite. 

Let then Vk = J- k [vo] denote the fc-fold application of 
the mapping J 1 " on a given initial potential v$ and choose 
y/a > C sup D. Then Eq. §TZ$) with a = C sup D/^/a im- 
plies llf/c+i — Vk\\ a < a fc ||^i — vo\\ a which means that 
the i>k are a Cauchy series. Since the set of potentials 
is a Banach space with the norm ||.|| Q [28[ and therefore 
complete, this series converges to a unique v , i.e. Vk — > v 
for k — > oo. According to our assumption the response 
function of Eq. (|15[) exists and hence q is functionally dif- 
ferentiate and consequently continuous as a functional 




FIG. 5: (color online). A simple example illustrating the 
meaning of the operator norm C[i>o, «i] which corresponds to 
the steepest tangent on the interval [«o,Wi]. 



of v. Therefore limj : _ ) . 00 q[vk] — q[v] which means that v 
solves Eq. (J5|) and hence is a fixed point. This establishes 
the existence of a Kohn-Sham system corresponding to 
the density n in Eq. (JB]) provided there is a supremum 
sup„ Cft'o,^ 7 ^]] when we range over potentials vo in a 
non-interacting system. 



VI. PERIODIC DENSITIES 

Let us summarize what we have found so far. The iter- 
ation Vk = T k [vo] will converge to an external potential 
giving a certain density n(xt) in an appropriately chosen 
a-norm. The mapping J- from the set of potentials 93 
onto itself depends on the chosen density via the solution 
of Eq. (fTU)) . Thus we have to invert the Sturm-Liouville 
operator in every iterative step. Only if the problem is 
regular at the boundary (see Appendix [B| we can invert 
in a straightforward manner by simply integrating twice. 
Otherwise we cannot right away fix the boundary condi- 
tions for nv' or v. For simplicity we now want to restrict 
our considerations in the following to the regular case. 

We will make the assumption of the regular case ex- 
plicit by assuming that the density is strictly positive, 
i.e. n{xt) > e > 0. This amounts to assume periodic 
boundary conditions for the quantum system. Then we 
can directly use Eq. (f25j) . It turns out to be convenient 
to define 

£(**)= / dy((yt)= I dy(q([v },yt)-d?n(yt)). 



(30) 
Then by partial integration and using periodicity we find 
from Eq. (gS]) that 



vi (xt) 



n(yt) 
where we defined the constant 



dy(d y K t (x,y))ayt) 
\y^ + c{t) [\y 1 



'(*) = 



dy 



1 



n(yt) 



1 r b 



dy 



n(yt) ' 

Z(vt) _ 

n{yt) ' 



(31) 



In order to eliminate the explicit dependence on q[vo\ in 
Eq. (f3"0")) we use the local force equation $5§ to write 

q([v },xt) = dln([vo\,xt) - d x (n([v ],xt)d x v (xt)). 

If we use this in Eq. (|3"Tj) we find 



v\ (xt) 



I -——ln([v },yt)dyv (yt) 



dzdt(n({vo], zt) — n(zt)) 



c(t) / dy 



i(yt) ' 



(32) 
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FIG. 6: (color online). The potential v(xt) in time from t = 
until T = 0.1. 



where 



m 



dy- 



i r b 



(yt)J J a n (yt) 



n{[vo\,yt)d y vo(yt) 



dzd?{n([v },zt)-n(zt)) 



Note that the form of Eq. ([52)1 is a very convenient one 
as it only involves densities and potentials, and there is 
therefore no need to calculate q([vo],xt) explicitly. Fur- 
thermore it is clear from the equation that the constant 
c(t) makes the potential periodic and that the explicit 
form fixes the gauge to v(at) = v(bt) = 0. 

Let us now give an explicit example of the iteration 
scheme. We take n(xt) to be the density produced in our 
previous example by potential of Eq. (J^H) (sec Fig. ©) 
with the initial state of Eq. (j2"5)l . Therefore our iteration 
scheme should recover potential (|29)l . We start the iter- 
ation with the initial guess vo(xt) = in the whole time 
interval. As the numerical inaccuracies tend to sum up 
in time we only look at a small grid and a short time 
interval. In principle we could perform the iteration for 
every time step and thus avoid the build up of inaccura- 
cies. However, here we arc not so much interested in a 
long-time propagation but in a proof of principle. Thus 
we start at i = and go only up to T = 0.1 atomic 
units. After one iteration we find (see Fig. (J7J) a first ap- 
proximation to the exact potential. The approximation 
becomes worse along the time axis. After 200 iterations 
we have almost converged to the exact potential in the 
first half of the time-interval, while in the second half nu- 
merical inaccuracies have build up (see fig. (J7J). That the 
iterative potentials rapidly approximate the exact poten- 
tial at earlier times is also obvious from the proof of the 
fixed-point, as the a- norm suppresses differences later in 
time strongly and thus the convergence is expected to 
be slower. As long as these numerical inaccuracies do 
not go out of hand, for instance by choosing the conver- 
gence time interval too large, then we find that after 1000 
iterations we have converged to a potential that is prac- 
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FIG. 7: (color online). The difference between the exact po- 
tential v(xt) and the iterated potentials vi(xt), V2oo(xt) and 
viooo(xt) in time from t — until T = 0.1. Note the change 
of scale. 



tically the same as v(xt) (see Fig. (J7J)). If we iterate even 
further we can make both potentials numerically indis- 
tinguishable. This little numerical illustration shows the 
convergence of the proposed iterative scheme and con- 
cludes the presentation of the fixed-point approach to 
the density-potential mapping in one spatial dimension. 
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VII. DISCUSSION AND CONCLUSION 

In this work we have given an extensive discussion of 
the density-potential mappings in quantum dynamics. 
We have derived in great detail all necessary equations to 
rewrite the density-potential mappings as a fixed-point 
question. We note, that a similar global approach was 
already introduced by Wijewardane and Ullrich in [29J . 
where they showed the numerical convergence of a time- 
dependent optimized effective potential calculation. The 
main equations are of Sturm-Liouville type and have been 
used already in the original work by Runge and Gross |7( . 
In order to uniquely solve these differential equations one 
needs to pose appropriate boundary conditions. These 
then fix the unique eigcnsolution to the eigenvalue zero, 
i.e. Eq. ([27]) . In the original Runge-Gross proof the po- 
tentials were assumed to vanish at infinity. If one allows 
for other boundary conditions then one can find differ- 
ent potentials leading to the same density, as has been 
shown in [30( . Therefore it is obvious that one can prove 
the Runge-Gross theorem only for a set of potentials with 
common boundary conditions. 

Further, we could extend the validity of the original 
fixed-point proof [15( ■ While the general linear response 
derivation in Sec. IIIII is independent of the dimensional- 
ity, we exploit the restriction to the one-dimensional case 
when solving the Sturm-Liouville boundary value prob- 
lems. In the case of periodic densities we can establish 
uniqueness and existence of a fixed-point for intcgrable 
potentials by direct integration. If the density becomes 
zero at the boundary we make use of singular Sturm- 
Liouville theory |16( . As shown in Appendix [D] if the 
density goes to zero at the boundary slower or equal to a 
quadratic function we can show uniqueness and existence 
of a fixed point for square-intcgrablc potentials. 

The discussed fixed-point approach to density- 
potential mappings provides a numerical scheme how to 
calculate the potential for a given density. In Sec. lVII wc 
show that the iterative sequence converges numerically to 
the exact external potential. We point out, that this pro- 
cedure has several potential applications. For instance, 
one can calculate the exact effective potential of Kohn- 
Sham TDDFT for a given interacting density. Hitherto 
this was only possible for special cases [31|. Further, one 
can use the density as the controlling functional variable 
in quantum control [32j. 

Finally one might be able to extend this fixed-point 
approach also to other functional theories, e.g . time- 
dependent current-density-functional theory |33( , lattice 
versions of density functional theory [34143 7j or supercon- 
ducting systems |35|. 
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Appendix A: Equivalence of the a-norms 

In this Appendix we will give a short summary of the 
functions spaces we are working in. We start by giving 
the mathematical precise form of the a-norm introduced 
in Sec. Hill 



|£ = ess sup 

te[t ,T\ 



*(*-*o) 



HtW 



Here the essential supremum is the supremum up to a 
set of Lebesgue-measure zero. For any a > this norm 
is equivalent to the norm 



II*' 



|g = csssup||w(t)|| p , 
te[t ,T] 



as can be seen by 

- a (T-t ) 



IHIS<IHIS<IHI- 



The space of functions which have finite ||-||o-norm is a 
Banach space, i.e. a complete normed vector space, de- 
noted by L°°([to,T], LP(I)) [H. Therefore the Banach 
space associated with any a-norm is isomorphic to this 
space. Consequently, if a sequence converges in some a- 
norm it converges in every a-norm. Thus all a-norms are 
equivalent and we can freely choose the constant a in our 
calculations. 

Further in the derivation of inequality (|13[) we used the 
operator norm C(i), i.e. Eq. (Tig)) . The variation therein 
(up to a normalization) goes over all g for which 



di||<7(i)|| p <oo. 



Those functions form the Banach space L p (Ho, T] , L P (I)) 
[281 ] . By the simple inequality 



'to 



dt\\g(t)\\ p < (T - t ) ess swp\\v 
te[t ,T] 



We can deduce that L°°([t ,T], LP (I)) c 

L p (jto,T],L p (I)) and thus we can include all func- 
tions with finite a-norm in our considerations. 



Appendix B: Sturm-Liouville boundary value 
problem 

In what follows we will give a brief sketch of self-adjoint 
Sturm-Liouville theory following the outline given in ref- 
erence [26| . For a thorough discussion of Sturm-Liouville 
theory we refer to [lq . 
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A general Sturm-Liouville boundary value problem 
reads as 



1 







w(x) \dx 



n{x) Tx 



+ k(x) <p{x) =\<p(x). (Bl) 



This is an eigenvalue equation for ip on an interval 
I = }a, b[ C R. In what follows we will use the so-called 
minimal coefficient conditions 

n _1 ,fe,ui G Ll oc (]a,b[), 

n(x) > and w(x) > on ]o, b[ a.e. 

Here "a.e." denotes "almost everywhere" , i.e. up to a set 
of Lebesque-measurc zero, and L\ oc {\a, b[) := {/ : / — >• 

C | ff dx\f(x)\ < oo V [a,P] C I}. Note that any 
continuous function is in Lj L oc (]a, b[). These conditions 
arc trivially fulfilled in our case, as we have w(x) = 1 
and k(x) = as well as a density n(x) which in general 
is continuous. In the following we will restrict ourselves 
to this special case. With this we will show, that we can 
always find a self-adjoint realization of the operator S = 
— d x \n(x)d x ] on the space of square-integrable functions 
L 2 (I) . We point out that an operator always consists of a 
"rule" , i.e. S, and a "domain" , i.e. which arc the functions 
it is allowed to act on. Actually, depending on I and 
n(x), we usually have an infinite number of self-adjoint 
realization which can be distinguished by different (self- 
adjoint) boundary conditions. We will introduce these 
boundary conditions in what follows. 

From the symmetry condition, i.e. (v\Su) — (Sv\u) = 0, 
the necessary self-adjoint boundary conditions have to 
guarantee that 

(v\Su) — (Sv\u) 
= lim {«(/?) [n(/3)V*(#>] " «*C8) HP)d p u{P)]} 

— lim {u(a) [n(a)d a v* (a)] — v*(a) [n(a)d a u(a)]} 



{u,v}(b)-{u,v}(a)=0 



(B2) 



for all v and u in the domain of the operator. Further 
the domain of any self-adjoint realization is a subset of 



D, 



{v g L 2 (I) \v, nv' e MiooQa, b[), Sv e L 2 (I) } , 



the so-called maximal domain. Here AC\ oc (]a, b[) is the 
set of locally absolutely continuous functions, i.e. v' exists 
a.e. and v(x) = v(a) + J dyv'(y) for all sub-intervals 
[a, 0\ C I and x £ [a, /?]. We note here, that the operator 
S becomes positive, i.e. (v\Sv) > 0, whenever 



v*(x) \n{x)d x v{x) 







for all v in its domain. This is a stronger restric- 
tion than sclf-adjointncss. Positivity is fulfilled, for in- 
stance, if one can choose homogeneous boundary condi- 
tions v(a) = v(b) = 0. However, one readily sees from 



the symmetry condition (|B2I) that the boundary condi- 
tions will in general not take such a simple form. In fact, 
the homogeneous boundary condition will in general only 
lead to a self-adjoint operator for the special case of a so- 
called regular Sturm-Liouville boundary value problem. 
In order to differ between the possible cases we introduce 
the following classification scheme: 

The lower endpoint a is called regular if 



a > 



and 



(l.r 
a n(x) 



< CO 



for an arbitrary c €. I. 
singular if cither 



The lower endpoint a is called 



a = — oo 



or 



dx 
a n(x) 



for an arbitrary c G /. If the endpoint is singular one 
cither has a limit- circle endpoint if for an arbitrary A G C 
and c G I any solution of Eq. (|B1[) obeys 



dx\^(x)\ 2 < oo, 



or one has a limit-point endpoint if for an arbitrary A G C 
and c G / at least one solution of Eq. (|B1[) obeys 



dx\ty(x)\ 2 = oo. 



Keep in mind that the equation is a second order dif- 
ferential equation and thus has two linearly independent 
solutions to every A G C. In a similar manner we have a 
classification scheme for the upper endpoint b. An exam- 
ple of two limit-circle endpoints is the well-known Leg- 
endre equation, i.e. n(x) = 1 — x 2 on / = ]— 1, 1[. Two 
linearly independent solutions to A = are 



and 



<P0,2( X ) 



VoA x ) = ! 



da; 1, fl + x 
n[x) 2 \1 — x 



With this it is easy to check the classification scheme. De- 
pending on the classification scheme of both, the lower 
endpoint a and the upper endpoint b we can pose differ- 
ent boundary conditions in order to have a self-adjoint 
operator S. 

The self-adjoint boundary conditions are directly re- 
lated to the condition (|B2[) and guarantee that the oper- 
ator S is symmetric. If the lower endpoint a is regular 
then one can pose the self-adjoint boundary condition 

Aiv(a) + A 2 n(a)v'(a) = 

with A 2 + A?, > 0. If the lower endpoint a is limit- 
circle then choose a pair /, g G D max with /,j€ K and 
{/, g}(a) 7^ 0. Subsequently 

A 1 {vj}(a) + A 2 {v,g}(a) = 
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with A\-\- A 2 , > is a self-adjoint boundary condition. 
These boundary conditions ensure that the term depend- 
ing on a in Eq. (|B2j) vanishes, i.e. {u, v}(a) = 0. We 
have according conditions for the upper cndpoint b. If 
the endpoint a is limit-point then no boundary condi- 
tion is needed nor allowed, as the normalizability of the 
functions becomes a necessary and sufficient condition 
to make the operator self-adjoint. Then for all func- 
tions in the maximal domain D max the boundary term 
{u, v}(a) = vanishes by construction. We note here, 
that in the case of regular or limit-circle boundary value 
problems we can also pose coupled boundary conditions, 
e.g. for a regular Sturm-Liouville boundary value prob- 
lem v(a) = v(b) and v'(a) = v'(b). For an example of 
limit-circle boundary conditions we again resort to the 
Lcgendre equation. A possible pair (/, g) is (vo.IjVo^) 
from above. If we choose A\ = 1 and A2 = for the lower 
cndpoint a = — 1 as well as the upper endpoint b = 1 we 
find the usual Lcgendre polynomials as eigenfunctions. 

It becomes evident that if we have two limit-point 
endpoints, there is only one possible self-adjoint domain 
for S, i.e. D max . Otherwise we have different possible 
choices for the domain of S. Irrespective of the choice of 
boundary conditions, our main interest lies in the spec- 
tral properties of the self-adjoint operator. Again we 
can rely on well-known facts from Sturm-Liouville the- 
ory. If both endpoints are either regular or limit-circle 
then we know that we have a pure point spectrum, i.e. 
only eigenvalues. And for separated boundary conditions 
we also know that we have simple eigenvalues, i.e. every 
eigenvalue has only one cigenfunction. In the case of 
coupled boundary conditions certain eigenvalues might 
have two eigenfunctions as can be seen from the nu- 
merical example in Sec. IIVI Depending on the bound- 
ary conditions it might occur that the eigenfunction to 
the eigenvalue zero is not the constant function. In or- 
der to assure tpo(x) = c we choose in the regular case 
either periodic, i.e. v(a) = v(b) and v'(a) = v'(b), or ho- 
mogeneous, i.e. v(a) = v(b) = 0, boundary conditions 
in accordance to the boundary conditions of our quan- 
tum system. For the limit-circle case we can always pose 
with the pair tpo^x) = c and ^0,2(2;) = / dy/n(y) the so 
called Friedrich's boundary condition A\ = 1 and A2 = 0, 



lim n(a)d a v{a) 



0, 



and accordingly for the upper endpoint. It is obvious 
that <po,2(x) cannot fulfill this condition while (fo t i(x) 
does. Therefore, for the regular as well as the limit-circle 
case the derivation in Sec. IIVI applies directly as we have 
a lowest non-zero eigenvalue and we can choose the zero- 
cigenfunction to be the constant function. 

If, however, one endpoint is limit-point (here automat- 
ically ipo(x) = c) one might have a continuous part in 
the spectrum. Then we have to make sure that the con- 
tinuous part is gapped away from zero. In the following 
Appendix [C] we will show that the continuum is indeed 
gapped away from zero in the case for the ground state 



density of a particle in a box by explicit calculation. Sub- 
sequently we will show in Appendix [D] that for a finite 
interval the continuum gap depends on the local behav- 
ior of the density near the boundaries. In particular, if 
the behavior of the density can be described as a power 
series with a lowest power p, the continuum is gapped 
away from zero for p < 2 and the gap closes for p > 2. 
How the inversion of the Sturm-Liouville operator with 
a spectral gap in the continuum can be defined is shown 
in Appendix [E] If the spectral gap closes we can not di- 
rectly apply the presented fixed-point approach on the 
set of square-integrable potentials. It seems reasonable 
to assume that in such a case the density can only be 
y-rcprcscntablc by a non-squarc-intcgrablc potential. 



Appendix C: Continuum gap for particle in the box 
ground state density 

In order to get some feeling for the continuum and the 
onset of the gap, we consider the lowest lying unnormal- 
ized density of one particle in a box, i.e. n(x) = cos 2 (x) 



on I = ]- 



7T 7T 

2 ' 2 



[. Then from 



dx ■ 



cos 2 (2) 



tan(x) 



00 



we deduce that a = — £ is a singular endpoint. An ac- 
cording calculation for the upper endpoint reveals that 
b = 7j is also singular. In order to check whether we 
have a limit-point or a limit-circle endpoint we need two 
linearly independent solutions to some eigenvalue A 6 C. 
Such two linearly independent solutions to the eigenvalue 
A = in this case are 






1. 



da; 

n(x) 



tan(x). 



(CI) 



We can then readily check from the classification scheme 
that due to 

/ dx tan 2 (2) = [tan(x) — x\ _ — > 00 

the lower endpoint a = — \ is limit-point. And from an 
according calculation we find that b = ? is a limit-point 
cndpoint as well. Therefore we find that D max constitutes 
the self-adjoint domain. 

Now we try to find the general solution for the differ- 
ential equation of the Sturm-Liouville problem 



d_ 
dx 



2, \ d 
cos {x) — 
ox 



<p(x) = Xip(x) 



First we try to eliminate the cos 2 (x) term by the following 
coordinate transformation y = cos 2 (x), so y <G [0, 1[. The 
differential equation then simplifies to 



' dy 2 



y( 1 ~y)i^-{ 2 y- 



d_ 
dy 



X 
Mi 



h(y) = 0, 
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where /i(cos 2 (x)) = ip(x). Note that we will only obtain 
the solution on half of the interval, since the inverse trans- 
formation is x = ± arccos y/y . The solution in the other 
half of the interval can be reconstructed by demanding 
continuity of the function and its derivative at x — 0. 

The new differential equation is almost the hyperge- 
ometric differential equation, except for the \/y term. 
This term can simply be dealt with by using the Frobe- 
nius method. We write h(jj) — y p f{y) and solve for 
p such that it annihilates the divergency, which gives 
pj = — i ± i\/l — 4A. The differential equation for f(y) 
reduces the following hypergcomctric differential equa- 
tion 



2/(1 - y)f"(y) + (2j>i 



2(pf + l)y)f'(y)- 
-p±(p± + l)/(y) 



with the solutions expressed in hypergcomctric functions 
2-F\(p A ,p x + 1) 2p A + §; y), so the full solution 
becomes 



ftiv) 



^( 



M 



(■>■) 



i F 1 (p±,p± + l,2p 



'■(x)). (C2) 



Note that for A < j the solutions <p^{x) are real and for 
A > t they are complex and the plus-minus solutions are 
each others complex conjugate, y A (x) = ip^(x). One 
can readily check that the general solution reduces to 
A = solutions found before (|C1[) . 

Now we found the solutions to the differential equa- 
tion, we turn back to the Sturm-Liouville problem. First 
we have to check which solutions to the differential equa- 
tion are normalizable to separate candidates for the point 
spectrum and the continuum spectrum. Since the hyper- 
geometric function is bounded over the interval, it is only 
the possible divergency of cos 2p > (x) at x = | in (|C2[) that 
can make the norm infinite. The small x 
the integrant is given as 



^ behavior of 



Iva(^)I s 



-1±V1-4A 



for A > \ 
for A < |, 



where we used that 2-Fi(o, b, c; 0) = 1. So only the func- 
tions <p~\( x ) with A < j are normalizable and could con- 
tribute to the point spectrum. 

To determine the point spectrum, we use that the func- 
tions should be smooth at x = 0. Since we effectively only 
solved the differential equation on half of the interval, say 
x G [0, -| [, this condition is not trivially satisfied. Since 
we have only one solution per eigenvalue, we can only 
construct a full solution using a ip~^ (x) with a G C for 
x < 0. Since we have two conditions to satisfy (continu- 
ity of the function itself and its derivative), we can only 
construct a solution if one of these conditions is satisfied 
automatically, irrespective of the value of a. Therefore, 
cither the value or the derivative needs to be zero at 



x = 0. The function values and derivatives of the general 
solutions at x = can be calculated to be 



^a(0) 



HiY(2p- 



M(0) = - 

ax 



Ap(p 



i) v^r(2 P + |) 



(C3) 



4p + 3 T(p + l)T(p + 2) 



Since T(2p + §) > for p > -±, we find that 0+(O) ^ 
for all A < j. However, for p — and p = — 1 the 
derivative vanishes at x — 0, so we find that </ 3 o"( a; ) is the 
only solution for A < j, such that the Sturm-Liouville 
operator is self-adjoint. Therefore, the only cigenfunction 
is the constant function with eigenvalue zero. 

Now we will determine which non-normalizablc solu- 
tions actually contribute to the continuum spectrum. Al- 
though they are not in D max , they are still required to 
give a self-adjoint operator. In particular, for / G D max 
the bracket {^a '/}(§) should vanish. The most diver- 
gent function in D max we can think of behaves near the 



boundary as x q , with q > 
we find 



Working out the bracket, 



WtJKi) 



{q-2pf) lim (/3-f 



n \2pf+q+l 

ivT=4A+ g +i 



( ff -2p±) lim_(/3-f 

(O for A > \ 
I oo for A < j, 



so we find that only functions with A > j build up the 
continuum spectrum. 

An alternative way to distill the continuum spectrum 
from the unnormalizable solutions comes from the spec- 
tral theorem. The spectral theorem states that the (gen- 
eralized) eigenfunctions of a self-adjoint operator form a 
basis for / G L 2 (] | , ^ [) , so in our case 



fix) = h + 



8=± 



d\f s (\)<p s x (x) 



thus the function / are considered as a wave-packet built 
from the continuum states (p\(x) and the constant func- 
tion. If we would have used the Laplace operator in ID 
with / = ]— oo, oo[, the continuum states would have been 
the plane waves, e ±ikx , and the integral would already 
start from and the expansion coefficients / would be 
the Fourier coefficients of /. Further note that the gen- 
eralized eigenfunction ip s x £ L 2 (] | , ^ [) , so they should be 
regarded as distributions. Therefore, the functions / S (A) 
have to be in the test-function space for the integral to 
be well defined. 

Since we required the generalized eigenfunctions to be 
such that the Sturm-Liouville operator is self-adjoint, 
the brackets ((B2]) between the generalized eigenfunctions 
should also vanish. In particular, for the upper endpoint, 
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b = t£, wc find in for A, A' > | 



2 ' 









0. 



This zero should be considered in a distributional sense, 
so if the bracket is integrated against test functions, the 
integral will vanish due to the infinite amount oscilla- 
tions near the edge, because the generalized cigcnfunc- 
tions near upper boundary behave as 



<M X ) 



1 



V 1 2 



-,± 5 \/4A^rL ln(x — § ) 



(C4) 



The vanishing of the integral can be formulated in 
a more precise manner by the Riemann-Lebesgue 
lemma. Similarly also for the lower endpoint we find 

It is now also rather obvious why the solutions ip x 
with A < 4 are not generalized cigenfunctions: they do 
not have the required infinite amount of oscillations to 
have a vanishing bracket. Indeed, if we check for the 
upper endpoint the symmetry condition (|B2[) for A > \ 
and A' < \ we find 



M^vHi) 






n .s2p*+2p-+l 



= 2(p A -p v ) km _(£-§) 

p— »• 2 
= 00, 

so indeed we recover that the solutions ip~^ (x) with A < j 
do not contribute to the continuum spectrum. 

The solutions for A > \ are complex. However, since 
also the complex conjugate has the same eigenvalue, we 
can combine them to construct real solutions. In par- 
ticular wc will choose the real solutions to transform as 
the irreducible representations of the symmetry of the 
problem, so separate them in gerade, (fx)'(0) = 0, and 
ungerade, y"(0) = 0, solutions. With the help of the 
values at x = (|C3|) . they readily constructed as 



^(x)=9f(^)'(0)S^+(x) 

-3%+)'(0) %>+(*), 
^(x)=sgn(x)(^+(0)%+(x) 

-Rp+(0) %>+(*)), 

where sgn(x) is the signum function. These generalized 
cigenfunction have been plotted in Fig. [5] for the eigen- 
values A = 1, 10, 100. When approaching the boundaries, 




FIG. 8: (color online). The non- normalized gerade (blue) and 
ungerade (red) generalized eigenfunctions for A = 1, 10, 100. 



the generalized eigenfunctions start to oscillate infinitely 
fast as mentioned before, due to the logarithmic term in 
the imaginary exponent (|C4[) . Therefore, the onset of 
the continuous spectrum is called the oscillation point, 
o"o = t- Thanks to the infinitely many oscillations and 
the divergence near the edge, the generalized eigenfunc- 
tions of the continuous spectrum are able to be orthogo- 
nal as required for symmetric operators, although there 
is an uncountable amount of them. In particular, we have 



(ipl\rt,)=Al6 ss ,S(\-\>), 



where s 



g,u 



and A\ is a normalization factor. The 



precise value of this normalization factor is important if 
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one wishes to compute the resolution of the identity 

1 f°° 1 

8(x-y) = -+J2 I d\—cp{(x)<p 8 x (y) 

~ -■ " T A 



s— g : u 4 



which is used in Appendix IE1 However, a spectral gap is 
all we need to validate our inversion and the associated 
inequalities. 



Appendix D: Continuum on a finite interval 

The example in the previous section shows that the 
continuum spectrum is actually a local property if we deal 
with a finite interval. The infinite amount of oscillations 
are pilled up near the edges where the density should de- 
cay sufficiently fast to zero such that the Sturm-Liouville 
operator does not blow up, so effectively supporting these 
oscillations. Assuming that the density decays as x p near 
the edges, we can indeed classify how fast the density has 
to decay to support a continuum and determine its onset. 

We will only focus on the lower boundary, since the 
results immediately carry over to the upper boundary. 
For convenience, we shift the interval such that the lower 
boundary is located at x = 0, so the interval under con- 
sideration will be ]0, c[, where c is some small positive 
number, such that n(x) = ax p is a good approximation 
to the real density. Without loss of generality, we assume 
that the boundary at c is regular and the boundary condi- 
tions could be smoothness of the solutions which depend 
on the precise form of n(x) over the whole interval. First 
we determine whether the boundary is singular 



d.v 
o n(x) 



dxx- p 



< oo for p < 1 
= oo for p > 1. 



Now we need to determine wether the singular cases are 
limit-circle or limit-point. For that we construct the two 
solutions at A = 



ipo,i(x) = 1 

<P0,2( X ) z 



d.v 
n(x) 




for p^ 1. 



The behaviour of the first solution is not problematic 
near the boundary for its normalizability. However, the 
second one might be problematic. In particular we have 



da; In (x) < oo 



,2-2p 



< oo for p < 



(l-p) 2 =oo forp>§. 



The results for the boundary classification have been 
compiled in Table |Hj 

To determine whether the singular limit-point cases 
support a continuum and where it starts, we need to 



TABLE II: End-point classification 



boundary classification 



-oo, 1[ 

[Lit 

[|,oo[ 



regular 

singular 

singular 



limit-circle 
limit-point 



solve the actual differential equation 

—d x [x p d x (p(x)] = \ip(x). 

The x p term between the derivatives can be eliminated 
by the following transformations 



V = hi(a;) 



y 



2\/A i 



-x 2 



(2-p) 



P 



for p = 2, 

for p 7^ 2. 



Since we have two different coordinate transformations, 
we need to deal with the cases p — 2 and p ^ 2 separately. 
Let us first consider the simplest one, p = 2. In this 
case the coordinate transformation turns the non-linear 
differential equation in a linear one which is straightfor- 
wardly solved by standard techniques. Transforming the 
solutions back, we find for the p = 2 case the following 
solutions 



2<pf(x) 



3 (-i±ivT=4A)ln(» 



r ±i v / T I: 4A 



Now we need to determine which solutions might con- 
tribute to the point spectrum. Checking the normaliza- 
tion gives 



dx\^(x)\ : 



< oo for ip+ <1/4 (x) 
= oo otherwise. 



Therefore, only the solutions ip^ <1 , 4 (x) can contribute 
to the point spectrum. The final selection depends on 
the boundary conditions at c and therefore, no more can 
be said about the point spectrum without additional in- 
formation. For A > \ we see that the solutions start to 



l 

4 

oscillate infinitely fast near the boundary, so we expect 
(To = \- Indeed, when we work out the bracket (|B2[) with 
x? 



and (7 < i we find 



{^,*«}(0)= lim( g - 



IzpiVl-4A).T 



g+iiiyT^iA 



ft) for(±,A>i) 
|oo for(-,A<i). 

So for a density decaying as n(x) = a^x 1 near the 
boundary, we find that there is a continuum starting at 
Co = «2/4. 

Now we turn to the cases p ^ 2. In these cases we 
are not so lucky that the coefficients in the differential 
equation become simply constants, but turns into 

h"(y) + -^-h'(y) + h(y) = 1 



py 
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where we defined h(y(x)^j = tp(x). Note the similarity 
with Bessel's differential equation; only the term in front 
of the first derivative is problematic. To eliminate this 
constant, we write the solution as h(y) = y a f(y) and 
choose a such that this constant becomes one. Following 
this strategy, one finds 

p-1 



and the equation for f(y) indeed reduces to Bessel's dif- 
ferential equation 

y 2 f'(y) + yf(y) + (y 2 -a 2 )f(y) = o. 

Performing all the back-substitutions, we can express the 
general solutions for p ^ 2 as 

^M = *»"-^| S |(E^-"> 

where J a (y) and Y a {y) are the Bessel functions of the 
first and second kind respectively. 

From these solutions we see that for p > 2, the Bessel 
functions start to oscillate infinitely fast when the ap- 
proach the boundary. Together with the divergence from 
the pre-factor they could constitute a continuum spec- 
trum. However, these oscillations are absent for p < 2, 
so we expect in these cases no continuum. Indeed, us- 
ing that for y ^> \a 2 — 1/4| the Bessel functions behave 
asymptotically as 



J a (y) ■ 
Y a (y) 



2 / 

cos [y 
iry 



sin I y 



~ ~ 4 
an n 
~2~ ~ 4 



we find that all solutions for p > 2 are not normalizable 
and that the brackets (|B2I) vanish for A > 0. In the case 
of p < 2 we need to use the approximation for small 
argument of the Bessel functions, < y <C \/a + 1, 



J<*(y) 

Y a (y) 



1 



r(o + i) 

'f(ln(y/2)+ 7 ) 



I» 



for a = 
for a > 0, 



where 7 ~ 0.5572 is the Eulcr-Mascheroni constant. Us- 
ing these asymptotic forms we find indeed that for p < | 
both solutions are square integrable and for | < p < 2 
only the p tp\ solutions are normalizable which can be 
used to construct the point-spectrum. Working out the 
brackets for the p tp\ solutions we find that they never 
satisfy the self-adjointncss condition (|B2[) . so there is no 
continuum spectrum if p < 2. 



TABLE III: Combined results of the lower boundary classi- 
fication and the onset of the continuum (oscillation point), 
o"o, assuming that the density decays as a p x p towards the 
boundary. 



p 


O"0 


point classification 


]-oo,l[ 


00 


regular point 


[l.i[ 


oc 


limit circle 


[|.2[ 


00 


limit point 


2 
]2,oo[ 


a p /A 



limit point 
limit point 



The results are summarized in Table IIIII Note that 
the situation with a finite continuum gap is rather ex- 
ceptional; it only occurs for p = 2. However, physically 
it is a very relevant one, since the density of a parti- 
cle in a box typically decays quadratically towards the 
boundary as we saw in Sec. [Cj Since the results are also 
valid for the upper boundary, we see that the side where 
the density decays fastest will determine whether there 
will be a continuum and the oscillation points, <tq. Fur- 
ther note that in the limit point - limit point case no 
additional boundary conditions are required/needed, so 
we can always make two linearly independent generalized 
cigenfunctions, i.e. the continuum is doubly degenerate. 
However, in the case of one limit point and one of the 
other boundary conditions we need to take a particular 
linear combination, so in that case the continuum will be 
simple, i.e. non-degenerate. 



Appendix E: Sturm-Liouville inversion with a 
continuous spectrum 



Here we perform the inversion of Sec. IIVI for a gen- 
eral self-adjoint operator, which might also have a con- 
tinuous part in its spectrum. We assume that there is 
a gap between the eigenvalue zero and the rest of the 
spectrum. From the spectral theorem for self-adjoint op- 
erators we know that every self-adjoint operator A has 
a unique spectral representation in terms of its spec- 
tral family (resolution of identity) E^ [lTj, which is an 
operator-valued function from R onto the set of orthog- 
onal projections. In physics one usually writes this res- 
olution of identity in terms of the Dirac notation, i.e. 
E A = J^du' |*^,)(*^,|. Here one can think of the 

|\fc£) as generalized eigenfunctions to the operator A. 
Thus we find 



A = 



AE*= duu\K)(*u\ 



If the spectrum only consists of eigenvalues, i.e. it is a 
pure point spectrum, then the integral becomes a sum 
over the discrete eigenvalues [l7| ■ The resolution of iden- 
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tity obeys 



as can be seen from 



d«|^)(^| = £, 



where i is the identity operator on the Hilbert space. 
In what follows we assume for notational simplicity that 
the Sturm-Liouville operator St = d x [n(xt)d x ] is posi- 
tive, i.e. its spectrum is in [0, oo[. The extension to the 
general case is straightforward. Then we can represent 
the inhomogencity of Eq. (fTCTj) as 

C(t)= / d«|¥ u >(¥ u |C(t)> 



du|* u )(*«IC(*)>+ / d«|* u )(*«IC(*)> 



where (,(xt) = g([vo] 5 xt) — d t n(xt) and \^ u ) are the gen- 
eralized eigenf unctions of St- Here we chose < e < Ai 
where Ai is a lower bound for the non-zero spectrum. By 
construction we know that £(£) is perpendicular to the 
A = cigenspace and thus we have 



§Mt)= I duu\* u )(* u \ I du'i|*„/)(^|C(i)) 



d«'-|^)(^|C(t))=C(*)- 



Here we used that {^f u \^u') = ${ u — u ')- Further we can 
deduce that v\(t) £ L 2 , since 



|| Vl (t)|| 2 = / d u -(C(t)|* M )(* M |C(t)) 

J e ** 

<I|| C (t)|| 2 <oo. 



In a similar manner we can then find 



m 



du|* u )(*«IC(*)> 



Therefore the solution to Eq. (fT0|) is 



ll«2(*)-«i(t)ll 2 <3llff(N.*)-?(N.*)ll 



«l(t) 



d« -!*„)(*„ |c(t)> 



If we then take D = 1/e we can derive inequality ([14 
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